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Abstract: Rapid industrialization and urbanization have led to the most serious habitat degradation in 
China, especially in the loess hilly area of the Yellow River Basin, where the ecological environment is 
relatively fragile. The contradiction between economic development and ecological environment 
protection has aroused widespread concern. In this study, we used the habitat quality of Integrated 
Valuation of Ecosystem Services and Tradeoffs (InVEST-HQ) model at different scales to evaluate the 
dynamic evolution characteristics of habitat quality in Lanzhou City, Gansu Province of China. The 
spatiotemporal variations of habitat quality were analyzed by spatial autocorrelation. A Geographical 
Detector (Geodetector) model was used to explore the driving factors that influencing the spatial 
differentiation of habitat quality, including natural factors, socio-economic factors, and ecological 
protection factors. The results showed that the habitat quality index of Lanzhou City decreased from 
0.4638 to 0.4548 during 2000—2018. The areas with reduced the habitat quality index were mainly located 
in the Yellow River Basin and Qinwangchuan Basin, where are the main urban areas and the new 
economic development areas, respectively. The spatial distribution of habitat quality presented a trend of 
high in the surrounding areas and low in the middle, and showed a significant positive spatial 
autocorrelation. With the increase of study scale, the spatial distribution of habitat quality changed from 
concentrated to dispersed. The spatial differentiation of habitat quality in the study area was the result of 
multiple factors. Among them, topographic relief and slope were the key factors. The synergistic 
enhancement among these driving factors intensified the spatial differentiation of habitat quality. The 
findings of this study can provide a scientific basis for land resources utilization and ecosystem restoration 
in the arid and semi-arid land. 
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1 Introduction 


The rapid industrialization and urbanization promote the continuous transformation and 
reconstruction of the spatial structure of land use, which has a great impact on the process of 
ecological landscape and ecosystem service level, and also threatens the security of regional 
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ecosystem (Zhou et al., 2014; Sallustio et al., 2017; Hao et al., 2017; Poniatowski et al., 2018; 
Kefalas et al., 2019). Specifically, land use change could affect the circulation of material flow 
and energy flow between habitat patches (Balkanlou et al., 2020; Shirmohammadi et al., 2020a, b), 
thereby changing the production capacity and service capacity of regional habitat (Aguilar et al., 
2019). Unreasonable land use patterns and rapid expansion of construction land had led to habitat 
degradation and fragmentation (Sallustio et al., 2017), which further reduced biodiversity (Krauss 
et al., 2010; Wilson et al., 2015). Therefore, strengthening the research on spatiotemporal 
variations and driving factors of habitat quality is of great value to ecosystem restoration and is an 
important way to construct the ecological security pattern, which has become a hot topic in the 
fields of ecology and geography (Sala et al., 2000; Newbold et al., 2015; He et al., 2017). 

Habitat quality is regarded as an important indicator of regional biodiversity and ecological 
service level (Maes et al., 2012; Moreira et al., 2018). It refers to the ability that the ecosystem 
provides suitable living conditions for the sustainable development of individuals and populations 
in a certain time and space. It is the foundation for the ecosystem to serve functions and provide 
services (Terrado et al., 2016; Zhu et al., 2020). According to the research object, data source, and 
research scale, there are mainly three methods to evaluate the evolution of habitat quality, such as 
traditional field survey of biodiversity and habitat (Vellend et al., 2008; Miller et al., 2009), 
ecological indicator assessment (Maes et al., 2012; Coates et al., 2016; Riedler and Lang, 2018), 
and ecological model (Costa et al., 2009; Terrado et al., 2016; Sallustio et al., 2017; Tang et al., 
2020; Zhu et al., 2020). Traditional field survey of biodiversity and habitat mostly focus on the 
spatial differentiation and influencing factors of single species and community habitat conditions 
in specific areas (Cardoso de Mendonga et al., 2003; Vellend et al., 2008; Miller et al., 2009). 
However, this method often need high data costs, and was only suitable for static analysis of 
habitat quality in a small region, and it is difficult to conduct research on the dynamic evolution 
of habitat quality (Vellend et al., 2008; Sun et al., 2019). Ecological indicator assessment is based 
on remote sensing to obtain ecological indices to evaluate biodiversity and habitat. Ecological 
indicators such as normalized difference vegetation index (NDVI) and net primary production 
(NPP) were analyzed frequently in recent study (Riedler and Lang, 2018). With the development 
of the Earth observation technologies, such as remote sensing (RS), geographical information 
system (GIS), and global positioning system (GPS), the widespread use of remote sensing image 
data, and the innovation of research methods, it is possible to quantify, evaluate, and simulate 
habitat quality at a large scale by using ecological models. Based on the habitat suitability index 
model, social value for ecosystem services model, and Integrated Valuation of Ecosystem 
Services and Tradeoffs (InVEST) model, researchers dynamically evaluated habitat quality at 
multiple scales (Zhu et al., 2020). Among them, the InVEST model is a mature and powerful tool 
for assessing habitat quality, because it comprehensively considers the suitability of habitat and 
anthropogenic threats to biodiversity, and provides more detailed information about the status of 
biodiversity (Terrado et al., 2016; Sallustio et al., 2017). In addition, the model is simple to 
operate, easy to obtain data, and has a strong ability to express results (Baral et al., 2014; He et al., 
2017; Sharp et al., 2020;). At present, the INVEST model has been widely used in habitat quality 
assessment at various scales, including national, provincial, transition zone, watershed, and nature 
reserves (Sallustio et al., 2017; Sun et al., 2018; Tang et al., 2020; Zhu et al., 2020; He et al., 
2022). 

Although there are many studies on the spatiotemporal evolution characteristics, influencing 
factors, and protection measures of habitat quality, there is still room for further research and 
deepening. First, most existing studies on habitat quality focus on the macro level or are limited to 
a certain scale, and there are few systematic and multi-scale studies on spatiotemporal evolution 
characteristics, driving mechanisms, and protection measures of habitat quality. In this paper, the 
"scale effect" was incorporated into the research process, breaking through the scale constraints of 
previous studies and revealing the transmission mechanism between different scales. Second, the 
research areas are mostly nature reserves, watersheds, provinces, countries, etc., but there are few 
studies on ecologically fragile and ecologically sensitive areas such as the loess hilly area in 
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natural transition zones. Third, a large number of researchers have focused on the spatiotemporal 
characteristics of habitat quality change and conservation measures, and there are few studies on 
driving factors of habitat quality. Most of the existing stidies on driving factors of habitat quality 
incorporate natural environment, geographical location, human interference, social economy, and 
other factors affecting habitat status into the single factor analysis model, such as scholars have 
analyzed the impact of urbanization (Tang et al., 2016; Li et al., 2018; McDonald et al., 2018; 
Song et al., 2020), land use changes (McDonald et al., 2009; Xie and Ng, 2013; Li et al., 2018), 
and human activities on habitat quality (Li et al., 2018; Sun et al., 2019), but there are few studies 
on comprehensively exploring the impact mechanism of habitat quality. In addition, most studies 
analyzed the spatiotemporal evolution mechanism of habitat quality based on the ordinary least 
squares (OLS) and geographically weighted regression (GWR) models (Li et al., 2017; Sun et al., 
2019), and few studies used the Geographical Detector (Geodetector) model to reveal the 
influence mechanism of the spatial differentiation of habitat quality. 

Lanzhou City, the capital of Gansu Province of China, is located in the loess hilly area of the 
upper reaches of the Yellow River, and the transitional zone of the monsoon region, arid and 
semi-arid region, and Qinghai-Tibet Plateau region, which is an ecologically sensitive and fragile 
region. Meanwhile, with the acceleration of social and economic development and_ the 
improvement of urbanization level in Lanzhou City, the fragile ecological environment and rapid 
urbanization process have led to the degradation of habitat quality, and have caused the decline of 
ecosystem services, as a result, the contradiction between socio-economic development and 
ecological protection has become more prominent. Thus, in this paper, we take Lanzhou City as 
the study area to fill the gap of habitat quality assessment in rapidly urbanization areas under the 
constraints of ecological fragility. It is expected to have practical significance for regional 
biodiversity protection and ecosystem restoration, and might provide a scientific basis for 
socio-economic development and ecological environmental protection. Specifically, the objectives 
of this study are to (1) apply the habitat quality of INVEST (InVEST-HQ) model to assess the 
evolution characteristics of habitat quality from 2000 to 2018; (2) analyze the spatiotemporal 
patterns of habitat quality; and (3) reveal the driving factors of the spatial differentiation of 
habitat quality. 


2 Materials and methods 


2.1 Study area 


Lanzhou City (35°34'-37°07'N, 102°35'-104°34'E) is located at the intersection of the Loess 
Plateau, Inner Mongolia Plateau, and Qinghai-Tibet Plateau of China, and is the geometric center of 
China's land area. The region is situated in a temperate continental climate zone with an average 
annual precipitation of 327 mm, an average annual evaporation of 1676 mm, an annual average 
temperature of 10.3°C, an average annual sunshine duration of 2447 h, and a frost-free period of 
180 d (Xu et al., 2021). The overall terrain is high in the northwest and low in the southeast, with an 
altitude range from 1418 to 3677 m. The landscape and vegetation types are diverse. The main land 
use types are grassland, cultivated land, and woodland. The Yellow River traverses the entire region 
from southwest to northeast, cutting through mountains and forming a beaded valley with canyons 
and basins alternately. The mountains from north to south face each other in the urban area, forming 
a belt-shaped valley basin with a length of about 35 km from east to west and a width of 2-8 km 
from north to south. There are three counties (Yuzhong, Yongdeng, and Gaolan counties) and five 
districts (Chengguan, Honggu, Xigu, Qilihe, and Anning districts) in Lanzhou City (Fig. 1). The 
total area of the city is 1.308x10* km?, and the permanent population of the city is 4.134x10° at the 
end of 2020 (Lanzhou Municipal Bureau of Statistics, 2021). 


2.2 Data sources and driving factors of habitat quality 


The land use data in 2000, 2010, and 2018 were downloaded from the Resource and Environment 
Science and Data Center of Chinese Academy of Sciences (http://www.resdc.cn/). The land use 
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Fig.1 Overview of Lanzhou City 


data was classified into 6 first-level categories, including cultivated land, forest land, grassland, 
water area, construction land, and unused land, and 24 second-level categories. 

The selection of driving factors of habitat quality is the key to the Geodetector model. In this 
study, we selected 12 driving factors based on comprehensive consideration of natural factors, 
socio-economic factors, and ecological protection factors that may affect the spatial distribution 
of habitat quality. The data sources of each driving factor are shown in Table 1. First, the data of 
each driving factor were resampled to a resolution of 100 mx100 m in ArcGIS software. Second, 
all the data were transformed into raster data and connected with the township (street) units in 
space. Finally, we classified and discretized the raster data according to the natural breakpoint 
method (Fig. 2). 


Table 1 Driving factors of habitat quality used in the Geographical Detector (Geodetector) model 
Primary factor Secondary factor Resolution/scale Data source 
Elevation (mm; X;) 30 mx30 m 


Natural factors 


Geospatial Data Cloud site, Computer Network 


Slop (°; X2) 30 mx30m Information Center, Chinese Academy of Sciences 
Topographic relief (m; X3) 30 mx30 m (http: wwwgseloud.cn) 
>10°C accumulated temperature 1000 mx1000 m 


(°C; Xa) 
Precipitation (mm; X5) 


Soil erosion (t/(km?-a); X6) 


Resource and Environment Science and Data Center 
1000 mx1000 m of Chinese Academy of Sciences 
(http://www.resdc.cn) 
1000 mx1000 m 


Socio-economic 
factor 


Population (X7) 


Gross domestic product (GDP) 
(x10° CNY; Xs) 


Lanzhou Statistical Yearbook (Lanzhou Municipal 
Bureau of Statistics, 2019) 


Road density (km/km?; Xo) 


1:1,000,000 National Geomatics Center of China (http://ngcc.cn) 


Distance from county (district) 
(km; Xio) 


30 mx30m Baidu map (https://map.baidu.com/) 


Nighttime light index 
(nW/(cm*+sr); X11) 


Resource and Environment Science and Data Center 


1000 mx 1000 (http://www.resdc.cn) 


Ecological 
protection factor 


Increased area of woodland (km?; 
X2) 


Resource and Environment Science and Data Center 
(http://www.resdc.cn) 


Note: — represents the secondary factor does not involve resolution or scale. 
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= 1706-1875 = 8.43-11.84 = 299-626 
= 1876-2015 = 11.85-14.62 = 627-835 
= 2016 2185 = 14.63-18.30 = 836-1173 
= 18.31-24.95 

= 2186-2480 = = 1174-1722 
(d) (e) (f) 


Soil erosion 


>10°C accumulated Precipitation (mm) (t/(km?sa)) 
temperature (°C) = 2315-2651 = 11.00-11.32 
= 13.07-22.14 = 2652-2959 = 11.33-11.90 
= 22.15-44.80 = 2960-3161 = 11.91-12.18 
= 44.81-59.78 = 3162-3287 = 12.19-12.39 
= 59.79-68.00 = 3288-3502 = 12.40-12.75 
= 68.01-83.00 = 3503-4207 = 12.76-13.48 
= 83.01-95.00 
(g) (h) (i) 
Population s Road density 
= 2904-10,627 GDP («10° CNY) (km/km?) 
= 10,628-19,579 oe = <0.15 
= 19,580-29,478 = 1.88-4.64 = 0.16-0.40 
= 29,479-45,225 = 465-837 may 
= 45,226-67,563 = 8.38-32.11 mol sae 
= 67,564-114,988 = >32.11 = >5.46 
(0) (k) (O) 
‘ Nighttime 
Distance from light index Increased area 
county (district) (km) (nW/(cm2.sr)) of woodland (km?) 
= <4,96 = 0.00-2.57 = —1.96—-1.28 
= 4,97-10.26 = 2.58-5.98 = ie ao 
= 10.27-16.80 = 5.99-12.86 aes 
= 16.81-27.06 = 12.87-25.23 = -0.06-0.79 9 20km 
= 27.07—40.18 = 25.24-45.71 = 0.80-5.91 
m >40.19 = 45.72-63.00 = 5.92-13.65 


Fig. 2 Classification of driving factors for the spatial differentiation of habitat quality in Lanzhou City. (a), 
elevation; (b), slop; (c), topographic relief; (d), =>10°C accumulated temperature; (e), precipitation; (f), soil 
erosion; (g), population; (h), gross domestic product (GDP); (i), road density; (j), distance from county (district); 
(k), nighttime light index; (1), increased area of woodland. 

2.3 Methods 

2.3.1 Habitat quality evaluation 

The InVEST-HQ model was used to assess the potential of ecosystems to provide survival and 
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reproduction for species (Sharp et al., 2020). Specifically, in this model, we set threat factors and 
habitat, and calculated the negative impact of threat factors on the habitat to obtain the 
degradation degree of the habitat, and then calculated the habitat quality through the suitability 
and degradation degree of the habitat. The habitat quality index is used to reflect habitat quality; 
as the index increases, the habitat quality and biodiversity increase (Li et al., 2018). The habitat 
quality is calculated as follows (Li et al., 2018): 


À 


D<. 
Q; =H ,|1-| —~ |}, (1) 
Dy +k“ 


where Qj is the habitat quality of raster x in land use type j; Hj is the habitat suitability of land 
use type j; Dx is the threat level of raster x in land use type j; z is a scaling parameter that reflects 
the spatial heterogeneity; and k is half the saturation constant, which is half of the maximum 
value of Dy. 


R wy, R 
Dy = pay eu / 2 Wingy Be jr > 2) 
day 
ig 1-[ » ) 6) 
: ee 
—0.299d 
is -e 7 = ; (4) 


where R is the number of threat factors; Y, is the number of grids of threat factor r in the land use 
map; y is the grid number of threat factors; w, is the weight of threat factor r (Table 2); r, is the 
stress value of grid y; i-xy represents the distance decay function, which can be expressed as a 
linear or exponential function of the distance from the threat factor to the habitat; Ø» is the 
accessible grid cell x, Sj, is the relative sensitivity of habitat type j to threat factor r (Table 3); dxy 
is the Euclidean distance between the habitat and the threat factor; dima, is the maximum 
interference radius of the threat factor r. 

Identifying threat factors to habitat is a key issue in the InVEST-HQ model. Based on relevant 
studies and the ecological characteristics of the study area (He et al., 2017; Li et al., 2018; Sun et 
al., 2019; Tang et al., 2020; Zhu et al., 2020), we selected urban land, rural settlement, railway, 
main road, cropland, bare land, and sandy land as the major threat factors to habitat quality (Table 
2). Besides, cropland (paddy field and dry land), forest land (forestland, shrub, sparse wood, and 
other wood land), grassland (high coverage grassland, moderate coverage grassland, and low 
coverage grassland), and waterbody (river and canal, lake, reservoir and pond, permanent glacier 
snow, mudflat, and beach) were selected as habitat types (Table 3). 

We calculated the maximum threat distance, the weight of each threat factor, and the relative 
sensitivity of habitat types to threat factors based on the InVEST User Guide (Sharp et al., 2020), 
expert scoring, and other relevant studies (Terrado et al., 2016; Tang et al., 2020). 


Table 2 Threat factors of habitat quality and calculated results from the habitat quality of Integrated Valuation 
of Ecosystem Services and Tradeoffs (InVEST-HQ) model 
Maximum threat 


Threat factor distance (km) Weight Decay type 


Urban land 6 1.0 Exponential decay 
Rural settlement 4 0.9 Exponential decay 
Railway 3 0.4 Linear decay 
Main road 3 0.6 Linear decay 
Cropland 1 0.3 Linear decay 
Bare land 2 0.2 Linear decay 
Sandy land 2 0.3 Linear decay 


Note: Decay type represents the type of impact of threat factors on habitat patches, including linear decay and exponential decay. 
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Table 3 Relative sensitivity of habitat types to threat factors 


Habitat type Urban land Rural settlement Railway Mainroad Cropland  Bareland Sandy land 
Paddy field 0.5 0.4 0.1 0.1 0.3 0.3 0.3 
Dry land 0.5 0.4 0.1 0.1 0.3 0.2 0.3 
Forestland 0.9 0.8 0.7 0.6 0.8 0.5 0.6 
Shrub 0.6 0.5 0.2 0.2 0.4 0.4 0.5 
Sparse wood 0.9 0.8 0.8 0.7 0.8 0.4 0.5 
Other wood land 0.9 0.9 0.8 0.6 0.7 0.4 0.5 
High coverage grassland 0.6 0.5 0.2 0.2 0.4 0.5 0.5 
Moderate coverage grassland 0.7 0.5 0.3 0.2 0.5 0.4 0.4 
Low coverage grassland 0.7 0.6 0.3 0.2 0.5 0.4 0.4 
River and canal 0.9 0.7 0.5 0.4 0.7 0.5 0.6 
Lake 0.9 0.8 0.6 0.5 0.7 0.6 0.6 
Reservoir and pond 0.9 0.8 0.6 0.5 0.7 0.5 0.6 
Permanent glacier snow 0.9 0.8 0.6 0.5 0.7 0.6 0.6 
Mudflat 0.8 0.7 0.5 0.5 0.6 0.3 0.4 
Beach 0.9 0.8 0.4 0.4 0.6 0.3 0.4 


2.3.2 Spatial autocorrelation analysis 


In this study, the global spatial autocorrelation index (Moran's J) was used to determine whether 
there is statistical agglomeration or dispersion in the spatial distribution of habitat quality. The 
calculation formula is expressed as follows (Zhu et al., 2020): 


n>, sae) 
Da haai -Yy 


where Moran's J is the global spatial autocorrelation index with ranging from —1.000 to 1.000; n is 
the total number of spatial units; Y; and Y; are the value of habitat quality of unit 7 and unit j, 


Moran's J = 


(5) 


respectively; and Y is the average value of Y;. 
2.3.3 Geographical Detector (Geodetector) 
In this paper, we used the factor detector (q value) to measure the interpretation degree of the 
independent variables to the spatial differentiation of habitat quality, and utilized interaction 
detector to reflect the interaction mechanism of the independent variables. The q value was 
calculated by Equation 6: 
1 L 

i are Ne ; (6) 
where q is the factor detector (q value); N represents the total number of spatial units; N; denotes 
the number of units in stratum i; and o? and o; denote habitat quality variance in population and 
stratum i, respectively. 


3 Results 


3.1 Spatiotemporal variations of habitat quality 


3.1.1 Variations of habitat quality at grid scale 

We used the InVEST-HQ model to obtain the spatial distribution pattern of habitat quality in 
Lanzhou City in 2000, 2010, and 2018 at grid scale (30 mx30 m), and divided the habitat quality 
index into five levels in ArcGIS, including level I (low), level II (relatively low), level II 
(medium), level IV (relatively high), and level V (high). 
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Table 4 showed that the average habitat quality index of Lanzhou City in 2000, 2010, and 2018 
was 0.4638, 0.4630, and 0.4548, respectively. The habitat quality declined slightly by 1.98% 
during 2000-2018. The proportions of habitat quality showed that the proportions of level I and 
level V increased gradually, while the proportions of level I, level III, and level IV presented a 
decreased trend, indicating that the polarization of habitat quality was obvious. 


Table 4 Changes of the habitat quality index in Lanzhou City in 2000, 2010, and 2018 


2000 2010 2018 
Level Range of habitat : : - 
ve quality index Proportion Average habitat Proportion Average habitat Proportion Average habitat 
(%) quality index (%) quality index (%) quality index 

I 0.0000—0.2000 2.98 3.79 6.10 

I 0.2000—-0.4000 57.28 56.48 54.39 

Il 0.4000—-0.6000 27.21 0.4638 27.16 0.4630 26.93 0.4548 

IV 0.6000—-0.8000 8.72 8.61 8.33 

V 0.8000—1.0000 3.81 3.96 4.25 


Figure 3 showed that the distribution pattern of habitat quality had significant spatial difference. 
The spatial distribution of level I was concentrated in the built-up areas of the Yellow River, 
Zhuanglang River, Huangshui River Valley, Yuanchuan River Valley, Yuzhong Basin, and 
Qinwangchuan Basin. The land use types were mainly urban land, rural settlement, and other 
construction land, which were affected by intensive human activities, and the habitat quality was 
fragmented. The level II and level III were mainly distributed in the loess hilly area of the Yellow 
River, the mountains in the north of Yuzhong County, and the hilly area of the western 
Zhuanglang River. The land use types were mainly cultivated land and grassland. The agriculture 


(a) 2000 


(c) 2018 


Habitat quality index mm 0.0000-0.2000 tm 0.2000—-0.4000 mm 0.4000—0.6000 
== 0.6000-0.8000 mm 0.8000—1.0000 


(d) 2000-2010 (e) 2010-2018 (£ 2000-2018 


Changes of the Changes of the Changes of the 


habitat quality index habitat quality inde: habitat quality index 
0.999905 0.999981 0.999981 
0.999982 0.999916 0.999982 


Fig. 3 Spatial distribution (a—c) and changes of the habitat quality index (d—f) in Lanzhou City from 2000 to 
2018 at grid scale. 
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and animal husbandry activity was more frequent, the ecosystem structure was unitary, and the 
ecological environment was fragile. The distribution of level IV and level V was covered by 
forestland and grassland located in the northwest and south mountains. As a result of high 
elevation, less human disturbance, and high vegetation coverage, habitat quality was generally 
higher. Above all, the spatial distribution of habitat quality has a strong spatial correlation with 
land use types. 

In order to further explore the spatial differentiation characteristics of habitat quality, we carried 
out the difference analysis of habitat quality and obtained the change maps of habitat quality in 
different periods (Fig. 3d—f). The results showed that habitat quality in Lanzhou City was stable 
from 2000 to 2018, but the local variation was obvious. Specifically, since 2000, with the 
acceleration of urbanization, the land on both sides of the Yellow River has been occupied by 
urban expansion, and the habitat quality has declined significantly. The improvement of habitat 
quality was mainly concentrated in the mountainous on the north and south sides of the Yellow 
River Basin and the north of the Qinwangchuan Basin, which was mainly due to the project of 
"Returning Farmland to Forest". After 2010, the areas with reduced habitat quality were mainly 
distributed in the plateaus at the edge of the Yellow River Basin and the southern part of the 
Qinwangchuan Basin. This is mainly due to the construction of a national-level new district (1.e., 
Lanzhou New District). The urbanization process has been accelerated obviously, and the river 
valley basin can no longer meet the demand of urban development, under the guidance of urban 
planning, cities began to break away from river valleys. The areas with improved habitat quality 
were mainly located in the mountainous areas in the northwest and south of the study area, which 
was mainly due to the restriction of the terrain, less human activities, and better habitat protection. 


3.1.2 Variations of habitat quality at town (street) scale 


Figure 4 showed that the spatiotemporal variation of habitat quality was more significant at town 
(street) scale. The towns (streets) with level I (<0.2584) habitat quality mainly concentrated in the 
main urban areas, and the number of towns (streets) increased gradually (26, 32, and 34 in 2000, 
2010, and 2018, respectively); the towns (streets) with level II (0.2584—0.4256) habitat quality 
formed four spatial distribution areas, and the number of towns (streets) decreased first and then 
increased (27, 24, and 25 in 2000, 2010, and 2018, respectively); the towns (streets) with level HI 
(0.4257—0.4723) habitat quality had a wider spatial distribution range but the number of towns 
(streets) decreased (32, 31, and 26 in 2000, 2010, and 2018, respectively); the towns (streets) with 
level IV (0.4724—0.5838) habitat quality were scattered in space distribution, and the number of 
towns (streets) first decreased and then increased (23, 18, and 21 in 2000, 2010, and 2018, 
respectively); the number of towns (streets) with level V (>0.5838) habitat quality was relatively 
small (3, 6, and 5 in 2000, 2010, and 2018, respectively), and they were mainly distributed in the 
northwest and southeast of the study area. 

Figure 4d-f showed the changes of habitat quality in Lanzhou City from 2000 to 2018 on the 
town (street) scale. From 2000 to 2010, 75 towns (streets) presented a declined habitat quality, 
among which 8 towns (streets) had a declined habitat quality with greater than 0.0600. There were 
36 towns (streets) with increased habitat quality, among which 7 towns (streets) increased by 
more than 0.0100. From 2010 to 2018, 74 towns (streets) showed a declined habitat quality, 
among which 14 towns (streets) had a declined habitat quality with greater than 0.0600. There 
were 37 towns (streets) with improved habitat quality, of which 16 towns (streets) increased by 
more than 0.0100. 

Generally speaking, the number of towns (streets) with a declined habitat quality was more 
than that with an increased habitat quality from 2000 to 2018. Among them, the towns (streets) 
with a large declined habitat quality formed two major agglomeration areas with Lanzhou New 
District and Lanzhou urban area as the center, which was mainly caused by the continuous 
expansion of urban space and the occupation of ecological land by construction land. The towns 
(streets) with a large increased habitat quality were mainly concentrated in the southwest of the 
study area, mainly due to the conservation of ecological land. 
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Fig. 4 Spatial distribution (a-c) and changes of the habitat quality index (d—f) in Lanzhou City from 2000 to 
2018 at town (street) scale. The classification criteria of the habitat quality index at town (street) scale is based on 
the natural breakpoint method. 


3.1.3 Variations of habitat quality at county (district) scale 

Figure 5 showed thevariations of habitat quality in different periods at county (district) scale. 
Habitat quality of each county (district) in Lanzhou City was relatively low. Among them, habitat 
quality in Anning and Chengguan districts was poor, while habitat quality in Qilihe District, 
Yongdeng County, and Yuzhong County was higher. In terms of different periods, from 2000 to 
2010, except for Qilihe District and Gaolan County, habitat quality of other counties (districts) 
showed a declined trend, especially in Anning District, where habitat quality decreased by 16.7%. 
From 2010 to 2018, habitat quality of all counties (districts) showed a downward trend. The 
change rate of habitat quality was —11.75% and —14.18% in Anning and Chengguan districts, 
respectively. The decline in the remaining counties (districts) was small (Fig. 5). 
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Fig. 5 Variations of the habitat quality index at county (district) scale in Lanzhou City from 2000 to 2018. Bars 
mean standard errors. 
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3.2 Spatial autocorrelation analysis of habitat quality 


In order to reflect the scale difference of spatial autocorrelation of habitat quality from both grid 
unit size and neighborhood radius, we constructed different grid units, including 1, 2, 5, and 10 
km, and used 1-2 times the grid unit size as the neighborhood radius to analyzes the spatial 
autocorrelation characteristics of habitat quality in Lanzhou City (Table 5). 

The results showed that habitat quality in 2000, 2010, and 2018 showed a significant positive 
spatial autocorrelation on each scale. With the increase of grid unit and neighborhood radius, the 
degree of spatial autocorrelation decreased, indicating that the spatial distribution of habitat 
quality was concentrated on the small scale and scattered on the large scale. In terms of the 
changes in spatial autocorrelation from 2000 to 2018, there were obvious differences on each 
scale. The spatial autocorrelation of the 1 km grid unit gradually decreased, the spatial 
autocorrelation of the 10 km grid unit gradually increased, and the spatial autocorrelation of the 2 
and 5 km grid units showed a trend of first decreased and then increased, indicating that the 
spatial distribution tends of habitat quality to be concentrated on the macro scale, more 
fragmented and discontinuous on the micro scale, and the meso scale has experienced a change 
process that tended to be concentrated after being dispersed. 


Table 5 Moran's J of habitat quality in Lanzhou City from 2000 to 2018 


Grid unit Neighborhood 2000 2010 2018 
(km) radius (km) Moran's! Z-value P-value Moran's] Z-value P-value Moran's] Z-value P-value 
i 1 0.774 125.286 0.000°° 0.779 126.208 0.000°" 0.797 129.005 0.000% 
2 0.734 184.491  0.000** 0.733 184.423 0.000%" 0.753 189.454 0.000% 
2 0.731 59.874 0.000°° 0.729 59.694  0.000** 0.745 60.989 0.000" 
‘ 4 0.687 87.465 0.000" 0.679 86.523 0.000°" 0.697 88.807 0.000" 
5 0.720 23.846 0.000°" 0.709 23.468 0.000" 0.715 23.670 0.000" 
i 10 0.629 33.010 0.000** 0.616 32.311 0.000°° 0.621 32.587 0.000%" 
iD 10 0.558 9.910  0.000** 0.548 9.740 0.000°" 0.542 9.606 0.000" 
20 0.395 11.430 0.000" 0.381 11.020 0.000°" 0.389 11.249 0.000" 


+ 


Note: “™* indicates significant at the confidence level of 0.01. 


3.3 Spatial differentiation mechanism of habitat quality using Geodetector model 


The above analyses found that there were obvious spatial differences in the evolution of habitat 
quality in Lanzhou City. Therefore, based on the township (street) scale, we selected 12 driving 
factors as independent variables from the aspects of natural factors, socio-economic factors, and 
ecological protection factors, and set the habitat quality index as dependent variable. Then, we 
used the Geodetector model to analyze the contribution rate and interactive detection results of 
driving factors, and to reveal the dominant driving factors of the spatial differentiation of habitat 
quality evolution in Lanzhou City and the interaction mechanism between the driving factors. 


3.3.1 Contribution rate of driving factors for the spatial differentiation of habitat quality 


Table 6 showed that the g value of the topographic relief (X3) and slope (X2) was 86.3% and 
76.8%, respectively, which are the dominant factors for the spatial differentiation of habitat 
quality. Secondly, the q value of the nighttime light index (X11) was 70.9 %, which means that 
human activities have a deep impact on habitat quality. Thirdly, the q value of soil erosion (X6), 
distance from county (district) (X10), gross domestic product (Xs), elevation (Xi), precipitation 
(X5), =10°C accumulated temperature (X4), population (X7), and road density (Xo) was between 
30.0%-50.0%, indicating that it has an important role in the spatial differentiation of habitat 
quality. In addition, the q value of increased area of woodland (X12) was only 14.5%. This 
indicates that increasing the ecological land area plays a certain role in improving habitat quality, 
but has a weak impact on the spatial differentiation of habitat quality. 
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Table 6 Contribution rate of driving factors for the spatial differentiation of habitat quality in Lanzhou City 
XxX, Xa Xa Xa Xs X6 X Xs Xo Xio Xu Xv 


q value 0.460 0.768 0.863 0.387 0.427 0.490 0.379 0.476 0.347 0.479 0.709 0.145 
P-value 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.087 


Note: q value is the factor detector, which indicates the interpretation degree of the independent variables (Xi—Xj2) to the spatial 
differentiation of habitat quality. 


3.3.2 Interaction effect of driving factors on the spatial distribution of habitat quality 


Table 7 showed the interactive detection results of driving factors for the spatial differentiation of 
habitat quality. The results indicated that the spatial differentiation of habitat quality in Lanzhou 
City was not caused by a single factor, but the result of a combination of multiple driving factors. 
The interactions of natural factors (Xi—X6) with natural factors and socio-economic factors were 
stronger than the interaction between socio-economic factors (X7—X11) and socio-economic factors. 
Among them, the interactions of topographic relief (X3) and slope (X2) with other factors were 
significantly stronger than the interaction between other factors. The results indicated that the 
synergistic enhancement of topographic relief, slope, and other factors had a great impact on the 
spatial differentiation of habitat quality in Lanzhou City. 


Table 7 Interactive detection of driving factors for the spatial differentiation of habitat quality in Lanzhou City 


Driving factor Xı X X3 Xs Xs X6 X; Xs Xo Xio Xu 
Xə 0.874" 
X3 0.886" 0.927" 
X4 0.532# 0.878# 0.882” 
Xs 0.650* 0.866" 0.893%  0.593* 
Xo 0.744" 0.841" 0.881" 0.747# 0.730" 
Xy 0.678" 0.847" 0.896" 0.660” 0.628" 0.668" 
Xs 0.644* 0.853# 0.906" 0.636* 0.651* 0.727*  0.548* 
Xo 0.592* 0.826" 0.886% 0.551* 0.555* 0.727* 0.621* 0.682" 
Xio 0.671# 0.870# 0.929# 0.667* 0.589% 0.710% 0.575# 0.638# 0.625% 
Xi 0.753" 0.898" 0.900" 0.755# 0.743* 0.771* 0.742# 0.740% 0.739* 0.738" 
X 0.579* 0.839% 0.886" 0.521% 0.548% 0.696* 0.571 0.621* 0.507** 0.642" 0.748" 


Note: * denotes interaction enhancement of two interaction enhancement and * denotes the nonlinear enhancement of two 
interaction enhancement. 


4 Discussion 


4.1 Spatiotemporal variation characteristics of habitat quality 


The results of this study showed that the habitat quality index in Lanzhou City decreased by 1.98% 
from 2000 to 2018. The main reason is that the area of construction land increases rapidly, 
resulting in the damage of habitat patches and the gradual decline of habitat quality. Habitat quality 
decreased mainly in the Yellow River Basin, Yuzhong Basin, and Qinwangchuan Basin. The spatial 
distribution of habitat quality showed a significant trend of high in the surrounding areas and low 
in the middle. The results are consistent with the study conducted by Xu et al. (2021) on the 
analysis of spatiotemporal variation characteristics of habitat quality in Lanzhou City. The overall 
level of habitat quality in Lanzhou City was low, and presented a continuous decline in time, with 
significant heterogeneity in space (Xu et al., 2021). In addition, the results in this study are similar 
to that in Liu et al. (2018). The areas with the low habitat quality index were mainly distributed in 
the regions with frequent human activities such as Yuzhong Basin, and the areas with the high 
habitat quality index were mainly distributed in higher altitude regions such as Xinglongshan 
Natural Conservation Area (Liu et al., 2018). This study also confirmed that habitat quality showed 
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significant positive autocorrelation at different spatial scales. With the increase of grid unit and 
neighborhood radius, the degree of spatial autocorrelation decreased. This indicated that the spatial 
distribution of habitat quality was concentrated on a small scale and scattered on a large scale. 


4.2 Complexity of driving factors for the spatial differentiation of habitat quality 


The spatial differentiation of habitat quality is the result of the interaction of multiple factors. In 
this study, the Geodetector model was used to analyze the influencing mechanism of driving 
factors on the spatial distribution of habitat quality in Lanzhou City. The results showed that 
natural factors determined the spatial distribution pattern of habitat quality, and socioeconomic 
factors and ecological protection factors were the main driving forces for habitat quality changes. 
These results about the influencing mechanism of habitat quality in Lanzhou City are consistent 
with those of Xu et al. (2021). The deterioration of habitat quality in the study area is mainly 
because of the increase of urban construction land in the process of urbanization. However, the 
improvement of habitat quality is mainly due to the construction of ecological function zones (Xu 
et al., 2021). In addition, the research on the mechanism of habitat quality should be further 
strengthened. Firstly, the relationship between drivers affecting habitat quality is relatively 
complex, it is difficult to obtain and quantify socio-economic data, and there are uncertainties in 
the results. Secondly, from the research results of natural factors, socio-economic factors, and 
ecological protection factors, the driving factors that cause the spatial differentiation of habitat 
quality are multiple and complex, but the coupling driving mechanism of multiple factors is still a 
problem. 


4.3 Potential habitat risk analysis in the study area 


Lanzhou City is regarded as a central city in northwest China. In recent years, due to the rapid 
development of industrialization and urbanization, significant changes have taken place in the 
spatial pattern of land use and ecosystem services (Lin et al., 2015; Hu and Pan, 2016; Yang et al., 
2021). Since 2010, with the construction of Lanzhou New District, the population in the area has 
been continuously agglomerated. In addition, unreasonable land use has reduced the areas of 
agricultural land and grassland, resulting in an increase in the degradation of regional habitat (Liu 
et al., 2016). These conclusions are consistent with the results of Zhou et al. (2017) using the 
ecosystem service value (ESV) method to evaluate the ecosystem service capacity in Lanzhou 
City. The total ESV in Lanzhou City decreased from 84.37x10° CNY in 1995 to 80.60x10* CNY 
in 2015, the most drastic changes were concentrated in the Yellow River Basin and 
Qinwangchuan Basin (Zhou et al., 2017). In addition, the downward trend of habitat quality in 
Lanzhou City has not changed, the potential habitat risk is relatively high, and the task of regional 
ecological protection is arduous. 


4.4 Suggestions and countermeasures for habitat quality protection 


In view of the status of habitat quality in Lanzhou City, we proposed the following protection 
measures. 

Lanzhou City is located in the loess hilly area of the Yellow River Basin, with a relatively 
fragile ecosystem. The forest area is small and the spatial distribution is uneven, mainly 
concentrated in the north and south of the study area, including important ecological function 
areas, such as Liancheng National Nature Reserve and Xinglongshan National Nature Reserve. 
Therefore, in the process of ecological restoration, we should increase the area of ecological land, 
and strengthen the protection of forbidden and restricted development zones, thereby improving 
ecological functions. Moreover, combined with the distribution of artificial corridors and the 
expansion characteristics of built-up areas, ecological corridors should be constructed between 
each core habitat patch to enhance the ecological effect of the corridor. Finally, we should 
strengthen the construction of the Yellow River Ecological Protection Belt, enhance the 
connectivity among Datong River, Huangshui River, Zhuanglang River, Wanchuan River, and 
other rivers, promote the transfer of material and energy as well as species exchange, diffusion, 
and migration between ecological sources in the study area. 
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Lanzhou City, as a typical river valley city, is restricted by topographic factors. Urban space 
and agricultural space are distributed along the river valley. The rapid expansion of construction 
land and the continuous occupation of cultivated land and habitat patches have led to the 
continuous fragmentation of ecological space, the gradual degradation of habitat quality, and 
continuous weakening of the connectivity of ecological corridors (Xu et al., 2021). Due to the 
rapid urbanization in the outer edge of the main urban area, Lanzhou New District, Yuzhong 
County, and other urban built-up areas, habitat quality has decreased significantly, the 
anti-interference ability of the ecosystem has weakened, and the contradiction between ecological 
protection and urbanization has intensified (Liu et al., 2016; Liu et al., 2018). Therefore, it is 
necessary to further strengthen the protection of ecological land, curb the rapid expansion of 
construction land, improve the efficiency of land use, maintain the dynamic balance between 
ecological protection and urban development, and promote the rational development of 
urbanization. 

In sum, through the assessment of habitat quality and the analysis of driving factors, it can 
provide reference for the optimization of land space pattern, scientific management and control of 
ecological space, and ecological construction and restoration. It is also helpful to alleviate the 
contradiction between urbanization and ecological protection in the upper reaches of the Yellow 
River Basin, and promote the coordinated development of economic growth and ecological 
environment protection. 


5 Conclusions 


Habitat quality is regarded as an important indicator of regional biodiversity and ecological 
service level. Evaluating the spatiotemporal evolution of habitat quality and its driving factors 
can provide a scientific basis for regional ecological protection and land use management. 
Based on the ArcGIS software, InVEST-HQ model, and Geodetector model, we analyzed the 
dynamic characteristics and spatial differentiation of habitat quality in Lanzhou City. The 
results showed that the habitat quality index in the study area decreased from 0.4638 to 0.4548 
during 2000-2018. The areas with decreased habitat quality were located in the Yellow River 
Basin and Qinwangchuan Basin, mainly due to the expansion of construction land. The spatial 
distribution of habitat quality showed an obvious trend of high in the surrounding areas and low 
in the middle, and showed significant positive spatial autocorrelation. As the grid unit and 
neighborhood radius increased, the degree of spatial autocorrelation decreased. This indicated 
that the spatial distribution of habitat quality was concentrated on small scales and scattered on 
large scales. The spatial differentiation of habitat quality in the study area was the result of 
combined effects of natural factors, socio-economic factors, and ecological protection factors. 
Among them, topographic relief and slope were the key driving factors of the spatial 
differentiation of habitat quality. Meanwhile, the interactions of natural factors with natural 
factors and socio-economic factors had a great impact on the spatial differentiation of habitat 
quality in Lanzhou City. 
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